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The ground state and its structure for a rotating, harmonically trapped iV-Boson system with a 
weak repulsive contact interaction are studied as the angular momentum L increases up to 3iV. 
We show that the ground state is generally a fragmented condensate due to angular momentum 
conservation. In response to an (arbitrarily weak) asymmetric perturbation of the trap, however, 
the fragmented ground state can be transformed into a single condensate state. We manifest this 
intrinsic instability by calculating the conditional probability distributions, which show patterns 
analogous to the boson density distributions predicted by mean-field theory. 
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Following the experimental realization of a dilute 
atomic Bose-Einstein condensate (BEC), the formation 
and properties of vortices in an atomic BEC have caused 
considerable interest both experimentally p|| and theo- 
retically (|-|| in the past few years. Although the re- 
cent demonstrations of vortex states by several different 
groups are in the Thomas-Fermi limit of strong inter- 
atomic interaction, a great deal of attention have also 
been attached to the nonvanishing angular momentum 
states of weakly interacting A^-Boson systems in har- 
monic traps. In Ref. ||, Wilkin et al. considered the 
case of an attractive interaction and showed that the 
ground state is uncondensed and is an example of the 
"fragmented^ condensate discussed by Noziere and Saint 
James @. Mottelson |j, Bertsch and Papenbrock ||] 
considered the lowest energy quantum states of a repul- 
sively interacting Bose gas when L ^ N. Wilkin et al. 
have employed a composite boson/fermion picture to de- 
scribe configurations beyond the one- vortex state S. A 
more tractable mean-field calculation performed by Butts 
and Rokhsar revealed a succession of transitions between 
stable vortex patterns of differing symmetries in the high 
angular momentum regime 0. The connection between 
the mean-field theory (MFT) and exact diagonalization 
scheme has been studied by Jackson et al. for a special 
case of L = 2N g. Finally, some analytical results have 
also been reported for the lowest energy states Q . 

In this Letter we address the question of whether the 
ground state of a weakly interacting A-Boson system 
with a given angular momentum is what one would nor- 
mally expect, i.e., a state with a single coherent Bose 
condensate, in which a mean-field approximation is valid. 
We propose that this is not the case and the ground 
state is generally a fragmented condensate in the pres- 
ence of the weakly repulsive interatomic interaction ex- 
cept L/N — or 1 in the thermodynamic limit. By eval- 
uating the macroscopic eigenvalues of the single-particle 
(SP) density matrix, we determine the degree of conden- 
sation. The origin of fragmentation turns out to be a 
requirement of the conservation of angular momentum. 



As a result, by turning on an (arbitrarily weak) asym- 
metric perturbation of the trap, the fragmented ground 
state can be easily deformed to a single condensate state 
PJ . This intrinsic instability can be manifested by the 
conditional probability distributions (CPDs) calculated 
for the ground state, which show patterns analogous to 
boson density distributions predicted by MFT. Note that 
the weakly interacting A^-Boson system considered here 
are quite similar to the spin-1 Bose gas studied by Ho 
and Yip [fl2f , in which the fragmentation originates from 
the spin conservation. 

We start from the model Hamiltonian describing N 
bosons in a two-dimensional harmonic trap interacting 
via a weak contact interaction. The SP spectrum is usu- 
ally expressed in terms of the angular momentum quan- 
tum number m and the radial quantum number n r , by 
En r . m = (2n r + \m\ + 1) fko. In the ground state of the 
system all the bosons are in states with n r = 0, and 
with m being zero or having the same sign as the total 
angular momentum. In the second quantized form, the 
Hamiltonian reads 

H = Ho + V, 

Ho = huj^2(j + l)aja,j, 
j 

V = \ 9 Vijkiafafa k ai, (1) 

i,j,k,l 

where Hq is the SP oscillator Hamiltonian and V is the 
two-body interaction between bosons. In the perturba- 
tive regime of weak interactions, Ng <C hu>. The opera- 
tor a, and af annihilate and create one boson in the SP 

j j 

oscillator state \ j) with energy {j + l)hw and angular mo- 
mentum jh, respectively, and obey the bosonic commuta- 
tion rules. The contact interaction elements are given by 
V m = S i+j , k+l 2-^{i + jy./iiljlkW) 1 / 2 g, and most 
of them are actually vanishing. For a given total angu- 
lar momentum L and number of bosons N, we consider 
the Fock space spanned by states \a) = \no,ni, ...,n/s) 
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with rij = N and ^jrij = L. Here rij denotes the 

j j 
occupation of the jth SP oscillator state \j) . There is a 
huge degeneracy corresponding to many different ways 
of distributing L quanta of angular momentum among 
N atoms. Here we restrict ourselves in a truncated Fock 
space of $C j $C j max = 6 flp^. To obtain the en- 
ergy spectra and the corresponding eigenstates, we set 
up the matrix elements in the Fock space basis, and sub- 
sequently diag onalize the matrix by using the Davidson 
algorithm fll4|| . 

A fragmented ground state. — First consider the SP 
density matrix in the form of 

p(r,r') = ^..^(r)p^(r') (2) 

% 3 

with ip m (r) = (r | m) . In Ref. (TJ] , Yang showed that the 
appearance of condensation is associated with the single 
macroscopic eigenvalue (i.e. of order N) of the density 
matrix p(r,r') with the "condensate wave function" be- 
ing the associated eigenvector, while the case of more 
than one macroscopic eigenvalue has been referred to as 
a "fragmented" condensate |H^]. The most important 
difference between the single and fragmented condensate 
is the lack of phase coherence of the latter. To find the 
eigenvalues of the SP density matrix, we write ]l5| ] 

Pij = Sp atpaf, (3) 

where the trace runs over all the N — 1 boson states, 
and the density matrix p = \^gs) {^gs\- It is readily 
seen that the eigenvalues are nothing but the occupation 
numbers of the SP oscillator state due to the conservation 
of the total angular momentum, namely p^ — SijHj. It is 
difficult to give a explicit expression for the occupation 
numbers rij. In the case of L = N, Wilkin et al. find 
that in the limit of iV -*■ oo, to the order 0(1 /N) Q, 

no = 1, n\ — N — 2, and n-i = 1. (4) 

They therefore conclude that the iV-Boson system is fully 
condensed into the one-vortex state in the thermody- 
namic limit. More detailed information can be obtained 
from the exact diagonalization calculations || . In figure 
1, we show the L dependence of the occupation numbers 
rij and their fluctuations Anj = (< n 2 > — < rij > 2 ) 1 / 2 
for j = 0, 1,2,3,4 for a system of TV = 40 bosons. When 
L ^ N the occupation numbers evolve rather smoothly 
as the angular momentum increases §|, while for L > N 
there are many kinks in the curves, reflecting the com- 
plexity of the ground states. The most prominent feature 
in the figure is that for a high angular momentum there 
are generally at least two significant occupation num- 
bers. For instance, at L = 70, the system has two large 
occupation numbers: rij ~ 9 and 23 for j = and 2, re- 
spectively. Evidently the case gives a "fragmented" con- 
densate. Although the present calculation is performed 
in the case of N = 40, the conclusion that a fragmented 
condensate ground state exists universally applies to the 



trapped, weakly interacting and rotating TV-Boson sys- 
tems with an arbitrary TV including the thermodynamic 
limit N -> oo @ . 

To examine the validity of the above statement, we in- 
vestigate the iV-dependence of the number of significantly 
occupied SP states by computing the inverse participa- 
tion ratio 0,11 : 

i c = Y.j^i N f- ( 5 ) 

The Iq is the first nontrivial moment of the distribu- 
tion of occupation numbers among the different SP states 
(note that Y2jrij/N = 1 by normalization). Its inverse 
1/Ic qualitatively measures the number of significantly 
occupied SP states. For example, Ic would be unity for 
a system that only had a single macroscopic occupied 
states; the maximum value of N is reached in the op- 
posite extreme, when all the SP oscillator states of the 
system are equally occupied. Figure 2 shows a plot of 
l/Ic as a function of angular momentum L for a sys- 
tem of N = 30,40,50 and 60 bosons. It is easy to see 
that in the regime of L/N < 1.6 the value of 1/Ic varies 
smoothly as L/N increases and shows little dependence 
on N. In particular, the variation of the peak height at 
L/N sa 1.6 is less than 3% as TV increases from 30 to 60 
(not shown in the figure). For L/N > 1.6, some irregu- 
lar small oscillations appear in the curves. However, the 
overall profile of 1/Ic is still nearly independent of N. 
These small oscillations are purely due to the finite N ef- 
fect JItJ and decay gradually with increasing N. One may 
expect them to vanish in the limit of TV — > oo. There- 
fore, we conclude that 1/Ic can be further used to quali- 
tatively measure the number of macro scopically occupied 
SP states in the thermodynamic limit, or in other words, 
to determine whether the ground state is fragmented or 
not. 

As shown in figure 2, there are two global minima (« 1) 
at L/N = and L/N — 1, which can be well inter- 
preted as a signature of single condensates. For other val- 
ues of L/N (especially in the high angular momentum), 
however, 1/Ic is generally larger than 2. This clearly 
indicates the fragmented nature of the corresponding 
ground states. Another notable feature in figure 2 is 
that the overall profile of 1/Ic exhibits a valley around 
L/N = 1.0, 1.8, and 2.4. This is consistent with the broad 
peaks of nj at L w 40, 70 and 90 for j — 1,2 and 3, re- 
spectively, as shown in figure 1. The number fluctuations 

1 /2 

Anj are in the order of 0(nj ) for these peaks, exhibit- 
ing a local characteristic of a single coherent condensate. 

The intrinsic spontaneous symmetry breaking of a frag- 
mented state. — Let us now consider the stability of such 
a fragmented ground state. In Ref. [ fulfill , the authors 
argued that the fragmented state is inherently unstable 
to the formation of a single condensate of well-defined 
phase. The essential idea is that even a weak perturba- 
tion that breaks the conservation laws will rapidly gen- 
erate phase coherence, modifying the density matrix de- 
terministically to give a unique macroscopically occupied 
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SP state. To support this point, we first show that the 
fragmented state and its corresponding single condensate 
state have the same energies in the limit of N — > oo, up 
to the order of O(gN). A similar conclusion has been re- 
ported by Jackson et al. for the special case of L = 2N 
H . In figure 3, the interaction energy Vint m units of gN 2 
is plotted as a function of L/N for a system of N — 20, 40, 
and 60 bosons. As N increases, Vint becomes closer in 
value to that of a single condensate V™/ , as predicted 
by the MFT @. The inset shows the energy difference 
&Vint — V™/ — Vi n t in units of gN. It is readily seen that 
all the AVi n t with different N are approximately located 
on a universal curve. This strongly suggests that AV, n t 
can be described by an approximate form 



AVnt = agN < fku, 



(6) 



where in the thermodynamic limit the factor a ~ 1 de- 
pends on L/N only and the inequality comes from our 
assumption of weak interaction. As a result, even a per- 
turbation of order 0(1 /N) can be enough to drive the 
fragmented state into a single condensate state. This fact 
clearly indicates that the fragmented state will sponta- 
neously break whatever the fragmentation was permitted 
by cylindrical symmetry in the first place. 

This result can be understood in another way by con- 
sidering the conditional probability distributions (CPDs) 
jUI that give the density correlation among bosons. We 
define the CPD for finding one boson at r given another 
at vq as 



V(v | v ) = 



*gs| Ei#j S ( r - r *)<5(vo - r,) |*gs> 



(N-i)(*Gs\'L i 8(yo-T J )\VGs) 



(7) 



Unlike the usual density distribution that is cylindrically 
symmetric under rotational invariant confinement, the 
CPD is asymmetric and reflects an intrinsic density dis- 
tribution 

What will an inherently unstable fragmented state 
evolve into if a weak perturbation is switched on? One 
may expect that the system will rapidly change into a 
state having the same intrinsic density distribution as the 
fragmented state, and simultaneously generate phase co- 
herence In view of this, the CPD gives the tendency 
of a system's evolution and can be regarded as a measure- 
ment of the possible spontaneous symmetry breaking. 

In figure 4, we show the L dependence of the CPDs 
for a system of N — 40 bosons. As expected, we ob- 
serve the successive vortex-like patterns of differing sym- 
metries, which are in good qualitative agreement with 
the mean-field calculations Q|. Both of them show a 
gradual transition for the formation of one- (fig. 4a) and 
two- vortex-like (fig. 4b) states in contrast to the rapid ap- 
pearance of the three- vortex-like state (fig. 4c) 0. As 
mentioned above, we identify this similarity as a signal 
of spontaneous symmetry breaking of fragmented states. 

On the other hand, one should not confuse CPDs with 
the "true" vortex patterns predicted by the MFT |t]] . The 



latter has phase coherence, which is not just well-defined 
in CPDs. Besides this, they have a different physical 
mechanism for the vortex emergence with the increasing 
L. For example, our results seem to show that the one- 
and two- vortex are produced at the center of the cloud of 
condensate, in apparent contradiction to the prediction 
of the MFT that the vortex enters the cloud from the 
low-density periphery. These differences may be resolved 
through the Josephson tunneling experiment suggested 
by Leggett and Sols |20|. Certainly, more accurate theo- 
retical studies on the fragmented state are required. 

In conclusion, we have studied the ground state of a 
weakly interacting .ZV-Boson system with a given angular 
momentum. We propose that the ground state is gener- 
ally a fragmented condensate state, which is rather frag- 
ile in response to a weak asymmetric perturbation. By 
calculating the corresponding CPDs, we manifest this in- 
trinsic instability. A comparison with the mean-field re- 
sults is also given. 
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Figures Captions 

Fig. 1. (color) Values of rij (in blue) and Arij (in red) 
of five lowest SP oscillator states as a function of L for 
a system of TV = 40 bosons with j max = 6. The cases of 
j = 5 and 6 are not shown due to their low occupancy. 

Fig. 2. 1/7 C versus L/N for a system of N = 30,40,50, 
and 60 bosons. Small oscillations at L/N > 1.6 are 
caused by the finite N effect. The overall profile of l/7<7 
is nearly independent of N, thus l/7c can be used to 
qualitatively measure the number of macroscopic occu- 
pation numbers in the limit of ./V — > oo. 

Fig. 3. (color) Vi n t in units of gN 2 as a function of 
L/N for a system of N — 20 (in green), 40 (in blue), and 
60 (in red) bosons. For comparison, the Vint predicted 
by MFT is also depicted by the dark solid line. The in- 
set shows the energy difference ISV int = V™J — Vint in 
units of gN. Note that all the AVmt with different N are 
approximately located on a universal curve. 



Fig. 4. Selected CPDs for a system of N = 40 bosons, 
(a), (b), (c) and (d) correspond to the emergence of 
vortex-like patterns with p-fold symmetry (p — 1,2,3, 4). 
In each panel, L increases in steps of one unit, and the 
starting value of L in (a), (b), (c), and (d) is 33,62,79, 
and 108, respectively. The values of x and y in each sub- 
plot range from —3.0 to +3.0. The given point vo is (0, 
1.0). For large N, the CPD is nearly independent of Vq. 
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